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Abstract 

>v 

The synthesis of helium in the early Universe depends on many input parame- 
^ ■ ters, including the value of the gravitational coupling during the period when the 

^ — , . nucleosynthesis takes place. We compute the primordial abundance of helium as 

function of the gravitational coupling, using a semi-analytical method, in order to 
track the influence of G in the primordial nucleosynthesis. To be specific, we con- 
■ struct a cosmological model with varying G, using the Brans-Dicke theory. The 

greater the value of G at nucleosynthesis period, the greater the abundance of he- 
' hum predicted. Using the observational data for the abundance of the primordial 

\f} • helium, constraints for the time variation of G are established. 
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The primordial nucleosynthesis is one of the most important achievements of the cosmo- 
logical standard model |l|, 0. Using the Einstein's equations, supposing a flat, isotropic 
and homogeneous space-time, considering a radiation dominated initial phase in the evo- 
lution of the Universe, it has been shown that it is possible, assuming initially an equal 
distribution of protons and neutrons, to obtain the primordial production of helium, lead- 
ing to an abundance of this element of about 24% of the mass of the Universe. The 
observational data indicates an abundance of helium of Y£ bs = 0.241 ± 0.0020 in mass 
fraction. Deuterium and lithium, as well as 3 He, are also produced. However, heavier 
elements can not be produced due to the absence of stable nuclei with A = 5 and 8; 
heavier nuclei are produced later in the stars. Taking the ratio of baryons to photons rj, 
as a free parameter, the predicted abundances agree with observations with a precision 
of some percents. For the moment, this is the earliest test of the standard model, giving 
confidence that the Universe followed, in general lines, the evolution predicted by the Hot 
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Big Bang model up to 1 second of age. It is possible that the analysis of the anisotropy 
of cosmic background radiation may lead to tests concerning earliest moments. But, this 
remains just a possibility, even if the results are becoming more and more consistents. 

The primordial nucleosynthesis, however, faces some controverses. As an example, 
it is considered now that lithium can also be produced in the stars, and the primordial 
abundance of this element could be smaller than it is assumed today ||. Moreover, 
recent measurements of the second peak of the spectrum of the anisotropy of the cosmic 
background radiation, made by the BOOMERANG and MAXIMA atmospheric balloons 
programs, indicate a density of baryons in the Universe that agrees only marginally with 
that necessary to have the correct primordial production of light elements. In fact the 
nucleosynthesis requires flbh 2 = 0.019±0.002 |5|], while the measurents of CMB anisotropy 
indicates flbh 2 = 0.032 ±0.005 ||. However, there are claims that this discrepancy is less 
important than indicated by those analysis [[7J. 

Even if the primordial nucleosynthesis remains one of the most impressive test of 
the standard model, it deserves yet more investigations. In particular, the primordial 
nucleosynthesis may be an arena to test the values of some input parameters. For example, 
the number of neutrino's families influence the primordial production of light elements. 
Hence, the primordial abundance of light nuclei permits to test if there are more than 
three families of leptons and quarks. 

In the present work, another effect on the primordial nucleosynthesis will be exploited. 
It will be verified how the value of the gravitational coupling affects the primordial pro- 
duction of light elements. Hence, the primordial nucleosynthesis will be taken as another 
way to verify if the gravitational coupling varies with time. In order to do this, we will 
work out a specific case: the Brans-Dicke theory, which is the prototype of a relativistic 
gravity theory with varying G, will be considered. The Brans-Dicke theory has an inter- 
esting connection with the low energy string action, this being another reason to study 
it. In the dust phase of the evolution of the Universe, power law type solution of the 
Brans-Dicke theory exhibits a decreasing gravitational coupling. This solution will be 
matched with the radiative solution, which in its simplest version coincides with that of 
the standard model. In this way the value of the Gr at the moment of the nucleosynthesis 
will be connected to its value today, G . The value of Gr depends on two input param- 
eters: the ratio of the densities of radiation and baryonic matter today; the Brans-Dicke 
coupling parameter uo. It will be shown that the greater the value of G at the moment of 
the nucleosynthesis, the greater the abundance of primordial helium. 

The precise calculation of the primordial abundance of light elements is a very hard 
task. It implies to use many numerical codes in order to evaluate the transmutation 
process involving protons and neutrons to obtain the final ratio between these nucleons, 
which determines the final abundance of helium. In order to have very precise predictions, 
radiative corrections to the rate of the reactions must be taken into account. In the present 
work we will adopt the semi-analytical method developed in ||. In this method the effects 
of the Fermi-Dirac statistic are neglected. This leads to a disagreement with respect to 
the precise calculation of the order of some percents. Since, the interest in this work it 
is to track the influence of a varying G in the calculation of the primordial abundance 
of light elements, it is interesting to sacrify somehow the precison in favour of analytical 
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expressions where the searched effects can be more easily tracked. 

This paper is organized as follows. In section 2, the Brans-Dicke cosmo logical model 
is developed, including the matching conditions. The relation between Gr and G is 
established. In section 3, the computation method presented in || is outlined, stressing 
the role played by G. In section 4, the primordial abundance for a varying G is obtained, 
in terms of the Brans-Dicke parameter u. In section 5 the conclusions are presented. 
Some considerations of how to reconcile a higher value for baryonic density, predicted by 
CMB anisotropy, with that one necessary to have abundance of primordial elements in 
agreement with observations, using a varying G model, are sketched. 



2 The cosmological scenario 

The Brans-Dicke theory incorporates the space-time variation of the gravitational coupling 
in a relativistic theory of gravity 0. In order to do so, it couples a scalar field non- 
minimally to gravity. Its lagrangian reads, 



-g 



4>R — uj—^ — 



+L m (1) 



where R is the Ricci scalar, is a scalar field connected with the gravitational coupling, 
uj is the Brans-Dicke parameter and L m is the Lagrangian of matter, which will be taken 
as a barotropic perfect fluid with an equation of state p = ap, — 1 < a < 1. General 
Relativity, and consequently a constant gravitational coupling, is recovered when uj — > oo. 
Using the flat Friedmann-Robertson- Walker metric, the equations of motion read, 

+ 3-0 = T7^V(P-3P) , (3) 
a 3 + 2uj 

p + 3-(p + p) = . (4) 
a 

There are two phases which interest us: the dust phase, a = 0, and the radiation phase 
a = 1/3. For these phases, the equations (@J|,|) admit the following power-law type 
solutions 



2+2uj 2 o 

— > a oc t 4 + 3 " , oc t 4 + 3 " , p m <x a ; (5) 
p = ^ — > a oc t 1//2 , 4> = constant , p r oc a -4 . (6) 

In the above expressions, p m and p r are the density for dust and radiation respect ivelly. 
For the radiative phase, the power-law solution coincides with that of the standard model, 
implying that the gravitational coupling is constant. But, since during the dust dominated 
phase, the gravitational coupling varies with time, the value of the gravitational coupling 
at the moment of the nucleosynthesis is, in this model, different from that normally used. 
Since the general features of the Brans-Dicke model for the radiative phase is the same 
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as in the standard model, we may compute the primordial abundance of light elements 
with a minimal modification. More general solutions, of course, can be found. We will 
discuss later the general features of these general solutions and their implications for 
nucleosynthesis. 

In order to extract precise predictions, the solutions for the dust and radiation phases, 
as well as their first derivatives, must be matched. The solutions are rewritten as 

a = a (t-t m )^ (p = 0) , a = b (t - t r ) 1/2 (p = |) , (7) 

PmO PrO / Q x 

Pm — —T- , Pr — —T , (P) 

a 6 a 4 

where ao, bo, t m , t r are constants; p m o and p r o are the density for matter and radiation 
today. Remembering that the gravitational coupling is connected with the value of the 

scalar field by G — f fijlfj j \ |10|| , and performing the matching of the solutions we find 

the following relation between the value of the gravitational coupling today Go and the 
value of the gravitational coupling during the radiative phase Gp. 

G R =(^)^G . (9) 

V PrO / 

It is more convenient to work with Planck's unity. Natural unities are employed: 
c = h — 1, G — 1 /y/Mp, Mp being the Planck's mass whose value today is Mp Q = 
1.221 x 10 19 GeV. Hence, 

Mp R J<^-)^Mp , (10) 
where Mpp is the value of the Planck's mass at the moment of the nucleosynthesis. 



3 The semi-analytical computation of the primordial 
nucleosynthesis 

The computation of the nucleosynthesis process in the early Universe involves three main 
steps. Initially, at very high energies, the ratio of neutrons to protons is equal to one. 
As the Universe expands, the temperature drops, and reactions involving the neutrons, 
protons, neutrinos, electrons convert neutrons into protons. At same time, since the 
neutrons are free, and unstable, they decay also into protons. This process continues 
until the energy is low enough in order the neutrons to be captured forming deuterium, 
from which the helium is formed. The quantity of helium synthetised in this process 
depends essentially on the quantity of neutrons that have survived up to the moment 
they are captured to form deuterium. The neutrons that are later used to form other 
elements like lithium are neglected in the present computation; anyway the others light 
elements besides helium represent a very small fraction of the total mass. The detailed 
analysis of all these process is quite involved, requiring the use of numerical codes to 
evaluate the different transmutation process. However, semi-analytical expressions can 
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be worked out if the Fermi-Dirac statistics is negletected. This has been done in ||, 
leading to values for the helium abundance with an error of some percents compared with 
the precise numerical calculation. 

Since all the evaluation of helium abundance in the approach used in || is very lengthy, 
we will just summarize the main steps and the relevant quantities. For details, the reader 
is invited to address himself to that work. First we define the ratio of neutrons with 
respect to the total baryon number: 



X(T) 



n n {T) 



n n (T) + n p (T) 



'ID 



where n n (T) and n p (T) are the numbers of neutrons and protons, respectivelly, as func- 
tions of the temperature T. The main equation controling this quantity is 



dX(t) 
dt 



X pn (t){l-X{t))-X np {t)X{t) 



(12) 



where X pn and X np are the rates of conversion of protons into neutrons and neutrons into 
protons respectively. The main process concerned are 



X np = X{y + n — > p + e ) + A(e + + n — * p + v) + X(n — > p + v + e ) 



(13) 



As an example, the first one is given by 



/■oo 

X(u + n~^p + e~) = A dp v p 2 u p e E e (l - f e )f v 

Jo 



(14) 



where A is a coupling constant, the p's denote the momenta of each particle involved in 
the process, E the energy and the /'s represents the Fermi-Dirac statistics factor. The last 
process in ([T3|) represents the neutron decay and it is not considered in a first evaluation. 
Later, the final results will be corrected taking it into account. 

In H the evaluation of the first two rates in (pT3|) is simplified by approaching the 
Fermi-Dirac statistics by the Maxwell-Boltzmann one. This is justified by the fact the 
temperatures concerned at the moment these process take place are smaller than the 
energies of the particles. After a lengthy evaluation of all process, we end up with the 
following expression for the neutron abundance factor: 



X(y) = X eq (y) + / dy'ey'xUy 1 ) ex V [K(y) - K{y')\ 



(15) 



with the following definitions: 

1 
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a = AAr{Am) 



Am 



(16) 
(17) 

(18) 
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where Am is the mass difference between neutrons and protons, Am = 1.294MeV. The 
fraction of neutrons to baryons at the end of all those process, X, is obtained by making 
V ->• oo, X = X(y -> oo). 

As stated before, initially the neutron decay is neglected. To correct the final results 
due to it we must evaluate the time taken in the capture process forming the deuterium. 
This time is given, in principle, by evaluating the capture process until a temperature of 
the order of the deuterium binding energy ~ 2.225MeV. However, lower energies (of 
order of E ~ O.lMeV) must be considered due to the fact that the enormous number of 
photons implies that deuterium dissociation continues to occur even when T 7 < eo- In 
peforming this analysis, we must take into account the evolution of the Universe, which 
in this case is reflected by the equation 



T 

i 

T-, 



8tt 



3 Ml 



(19) 



since during the radiative phase a oc 1/T 7 , T 7 being the photon temperature, which 
is approximately equal to the neutrino temperature at the relevant temperature scales 
considered in this computation. Moreover, the energy density is given by 



p = N ( 



7T 



eff 



rp4 

30 v 



where 



11\ 4/3 

N v +[—) iV 7 ~ 
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N eff = j 
The final results imply that the time of capture is given by 



tr. 



45 



167T 3 A^, 



eff 



1/2 



4 



11 



2 /3 Mr 



The constant t is 



11 



6N, 



eff 



:-) 1/3 - 1 

v 4 ; 



T 2 
1 70 



45 



16vr 3 iV e 



V2 /1U 2/3 Mp 



(20) 



(21) 



(22) 



(23) 



where T 7 o is the photon temperature when the neutron capture is accomplished. The final 
abundance is given by 

X f = exp (-t/r)X (24) 
and the final helium abundance by weight is 



Y d = 2X 



f 



(25) 



Applying all these steps to the standard model, it results Ya ~ 0.243. Note that this 
result is not exactly the same found in ||, which is Ya ~ 0.247. We attribute this 
small discrepancy to the fact that we used a Mathematica program instead of a pocket 
calculator. Moreover, it is point out in M that their results is a kind of upper limit using 
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the semi-analytical method described there. A more precise computation, using numerical 



code, gives Y 4 ~ 0.241 |3|. 



The above expressions resume briefly the main steps. More important, it has been 
shown explicitly where the value of the Planck's mass (consequently, the gravitational 
coupling value) appears. This will enable us to compute the variation of the helium 
abundance as function of Gp, and consequently as function of the Brans-Dicke parameter 

UJ. 

In H , it was computed the variation on the helium abundance due to variations of the 
chemical potential of the electron neutrino /i, the number of neutrinos N u , the neutron 
time-life r and ratio of photons to baryons rj. They find 

AY 4 = -0.25u + 0.0UAN v + 0.18— + 0.009 In — (26) 

where t]q = 5 x 10~ 10 is the adopted value for the ratio of baryons to photons. One of our 
goals is to add to this expression the term concerning the variation of the gravitational 
coupling as a perturbation around its value computed using the value of G today. 



4 Helium abundance with varying G 



Now we turn to the computation of the helium abundance with the cosmological model 
developed in section 2. First of all, we notice that since in the radiative phase the scalar 
field is constant, all development exhibits in the previous section is valid; we must just 
to compute the value of the gravitational coupling (Planck's mass) in view of the fact 



that in the later dust phase the gravitational coupling varies with time. We fix in ( [T0| ) 
Pmo ~ 10~ 29 g/cm 3 and p r o ~ 0.950 x 10~ 33 g/cm 3 , which are approximately the densities 
of matter and radiation today |TJj]. Using Mpo ~ 1.221 x 10 19 Gew, we can then compute 
the values of Mpp at the radiative phase in terms of to. With the value of Mpp, we can 
compute the helium mass fraction. The values are the following: 



UJ 


M PR (in GeV) 


Y 4 


-1 


oo 


1 


-0.5 


1.160 x 10 15 


0.961 





1.190 x 10 17 


0.818 


1 


1.205 x 10 18 


0.612 


10 


8.015 x 10 18 


0.318 


50 


1.115 x 10 19 


0.259 


100 


1.166 x 10 19 


0.251 


500 


1.210 x 10 19 


0.244 


1,000 


1.215 x 10 19 


0.243 


10,000 


1.220 x 10 19 


0.243 



From this table, we can see that all matter becomes composed of helium for uj = — 1 
(the case where the Brans-Dicke cosmology coincides with the string cosmology at low 
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energy level |15|]), and it approaches the asymptotic value of the standard model as u 
increases. This result is easy to understand. Since in the radiative phase the Fried- 
mann equation ( |I9| ) is valid, as G increases the rate of the expansion of the Universe 
also increases. Hence, the temperature drops very quickly and the duration of the nucle- 
osynthesis era becomes shorter. All effects described in the preceding section contribute 
to decrease the quantity of neutrons. If the duration of the nucleosynthesis era becomes 
smaller, than more neutrons survive to be captured in deuterium, forming helium later. 
When u) — — 1 the value of the gravitational coupling diverges, there is an instantaneous 
transition from the radiative phase to the dust phase, and all initial neutrons survive. 

We may evaluate now how a small change in the value of the gravitational coupling 
affects the nucleosynthesis as an approximation to the standard model. This will permit us 
to establish constraints on the variation of G by using the nucleosynthesis observational 
results. The value of G affects essentially the parameters b, t c and T 70 . Writing G = 
G + AG, (Mp = Mp + AMp), and introducing this quantity in the computation steps 
described before, we find that the helium matter fraction varies as 

AG 

AY 4 = 0.088— . (27) 
Go 

It is generally argued that the nucleosynthesis observational results coincides with the 
theoretical results by with a precision of 1 — 2%. Hence 

A G G i n 1 , 

„ i _> „ KT^years- 1 (28) 

where we have, in the last step, divided the first expression by the age of the Universe, 
assumed to be tu ~ 2 x 10 10 years, to obtain an estimation of the fractional time variation 
of G. This agrees in order of magnitude with the estimations of the time variation of 
G using other methods W\ O, |TE|. Of course, a more precise comparison with other 



experimental determinations of ^ is not possible due to the approximations made. 



5 Conclusions 

Primordial nucleosynthesis is considered as one of the most precise tests of the standard 
cosmo logical scenario. With it, the abundances of helium, deuterium and lythium, due 
to the primordial processes, are computed. Comparison with observations showed that 
the predicted abundances agree with the observed ones by some percents. This is an 
impressive result since the primordial nucleosynthesis ocurred in the first seconds of the 
existence of the Universe. This high precise cosmological test may be used to constrain 
the values of some input parameters, like the number of neutrinos and the ratio of baryons 
to photons. 

In the present work, we have estimated the influence of a possible variation of the 
gravitational coupling G in the predicted abundance of light elements. In order to be 
specific, a Brans-Dicke cosmological scenario was constructed, with a radiative phase 
identical to that of the standard model, followed by a dust phase. The gravitational 
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coupling decreases during the dust phase, while it remains constant during the radiative 
phase. Hence, the value of G at the moment of the nucleosynthesis is, in this model, 
higher than the usual one. 

The results indicate that for a larger G, the abundance of helium is higher. This 
is due to the fact that all the process occuring during the nucleosynthesis contribute to 
decrease the number of neutrons. A larger value of G leads to a faster expansion of 
the Universe, decreasing the duration of the nucleosynthesis era. Consequently, the final 
number of neutrons that will survive, forming later the helium, is higher. The precision of 
the measurements of the abundance of the primordial elements permits then to estimate 
the possible variation of G. We found however, that window of allowed values is not larger 
than that one obtained by another kind of experiment. 

In all computation, it was used the simplest power-law type solution of the Brans- 
Dicke theory for the radiative phase, which is identitical to the standard model one. This 
simplify a lot the evaluation of the helium abundance. However, it is possible in principle 
to consider more complicated scenarios. In fact, for the radiative phase, the Brans-Dicke 
scalar-field can be written as 

j>=- 3 , (29) 

where C is a constant. If C is positive, the gravitational coupling decreases with time 
during the radiative phase; if it is negative, the gravitational coupling increases with time. 
The case studied before corresponds to C = 0. For C ^ 0, we have the following solution 
for the scale factor: 



a(£) = a 



I 3 



I 3 



(30) 



where B 



§C,0 



lu + | and £ is the conformal time 



■ /* 

a 



(31) 



It will be interesting to consider such more general model, including a possible increase 
of the gravitational coupling during the radiative era, since it can be a possible solu- 
tion to the discrepancy of the baryonic density obtained from nucleosynthesis and from 
CMB anisotropy. This leads, however, to more important modifications in the preceding 
computation, since some basic equations, like (|19"|) must be modified by including the con- 
tribution of the non-minimal coupled scalar field. We hope to present this more general 
analysis in the future. 
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